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1. Introduction 

Integrable Quantum Field Theories in 2 dimensions are successfully studied by using the S- 
matrix approach both for massive [p and massless theories [|2|, On the other hand the study 
of non-integrable models is complicated because their scattering amplitudes are generally non- 
elastic. These models are in fact characterized by particle production, resonances, decay events, 
etc. along with the simple pole structure that features integrable theories. 

The 0(3) non-linear a model in 2 dimensions with a 8 term is a well-known example of 
non-integrable theory. It is defined by the action [||] 



Y J d 2 x(d^{x)f -id J d 2 xQ(x) , (1.1) 



where g is the coupling constant, <j>(x) is a 3-component unit vector and Q(x) is the topological 
charge density operator ^ (a,b,c are group indices that run from 1 to 3 and n, v are 2-dimensional 
spacetime indices) 

Q{x) = ^ v e ahc 4 a (x)d^\x)d v f(x) . (1.2) 

The integration of Q(x) throughout the whole 2-dimensional spacetime yields the total topological 
charge Q which takes on integer numbers that reveal the winding of configurations over the sphere 
S 2 . Topological charge (— )1 configurations are called (anti)instantons. 

When 8 = or % the model is integrable and its spectrum is as follows. For 8 = it contains a 
triplet of massive scalars whose mass has been analytically calculated [Q] and numerically verified 
within a 2%-3% error [0]. At 8 = n it has been conjectured that the theory becomes massless ® B 
(for a numerical analysis of the corresponding universality class see [§]). 

The study of the evolution of the spectrum as 8 moves from % towards lower values is worth- 
while. The model contains a triplet and a singlet states whose masses (mj and ms respectively) 
are proportional to (tc — 8) 2 ^ 3 for < (71 — 6) <C 1 dlO |. By using Form Factor Perturbation The- 



ory fliT] ] it has been shown that the singlet is heavier than the triplet [jl2|] (see also Ref. [jl3|]). In 
particular at 8 rj n the ratio of singlet to triplet masses is ms / mj \/3. 

On the other hand it is known that at 6 = there is no lingering trace of the singlet state. This 
fact suggests that the singlet mass may diverge as 6 approaches zero, thus decoupling from the 
whole theory. Consequently, by continuity arguments, it seems plausible that there exist a critical 
d c where the singlet mass becomes exactly twice the triplet mass in such a way that for all 6 < 9 C 
the singlet decays into states belonging to the triplet. 

In Fig. |] a sketch representing the qualitative behaviour expected for the 6 dependence of the 
spectrum of the theory is shown. We are currently making a Monte Carlo study of the model to 
examine some of the salient features shown in this Figure. In the present progress report we give 
clear numerical indications that mj{6) vanishes at the endpoint O en d = % (for an extensive account 
see [0]). 

To this end we have prepared a code that simulates the 2-dimensional 0(3) sigma model 
with an imaginary 8 term. After extracting the triplet mass with this code, one can analytically 
continuate the results to real values of 8 thus obtaining a numerical estimate of the lower curve 
in Fig. |l|. We focus our attention on the location of the endpoint 8 en d on this curve. In section 2 
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Figure 1: Sketch of the hypothesized 6 dependence for the singlet (dashed line) and triplet (continuous line) 
masses in the 0(3) non-linear sigma model in 2 dimensions. C is the value where the singlet mass becomes 
larger than twice the triplet mass, iuhmn = mj(Q = 0) is the mass calculated in Ref. [§]. e nd is the endpoint 
where the theory becomes massless. It has been conjectured that enc i = n. 

we shall shortly describe the updating algorithm used at imaginary 6. In section 3 the analytic 
continuation is performed and the main results are exhibited. 



2. The Lattice implementation 

We have regularized the model in Eq.( |l.l[ ) by using the standard lattice action 

S L =A L -i6 L Q L , A L = - j 8£0(x)-0(x + J u), Gt = £fii(jc), (2.1) 

X,fl X 

where /3 is the inverse bare lattice coupling constant for this standard regularization and 6l is the 
bare theta parameter. The corresponding endpoint and critical point will be called 0L.end and d L c . 



The lattice operator for the density of topological charge is defined as in []15[], 
1 



Ql(x) 



327T 



e^ v e abc r (x) U b {x + ju) - 0"(x - £) )U c (x + v) - c (x - v 



(2.2) 



The coordinate sites are labelled by x = (x\,X2)- The total topological charge Qi obtained from 
the sum over the whole lattice volume of the above expression does not yield integer numbers on 
single configurations. This is not an inconvenience inasmuch as Quantum Field Theory deals with 
quantum averages of renormalized operators over many configurations. Thus the above definition 
for the density of topological charge ( [2J| ) must be renormalized, the renormalized charge being Q 



Ql = Z q Q, (2.3) 
and Zq being the corresponding renormalization constant. It can be calculated either by perturba- 



tion theory [16] or by a non-perturbative numerical method [17]. The meaning of this constant 
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is simple: it accounts for the average over quantum fluctuations in such a way that Q yields an 
integer value. Actually this last observation is the basis for the above-mentioned non-perturbative 
numerical method (see below). 



The reason to choose the expression ( ]2.2| ) to be used in the lattice action is that, as we shall see 
later, it allows the introduction of a variant of the cluster algorithm in the presence of the 6 term. 
There exist other lattice regularizations of the operator Q(x) that do not require the computation of 
Zq (since it is exactly 1 for all j8); however it is not easy to introduce a fast cluster algorithm for 
them. 

Our interest concerns the spectrum of the theory for varying 6. This parameter is related to the 
corresponding bare 6l by the expression 6 =Zq 6l. Therefore the value of the bare theta parameter 
$L.end where the mass vanishes is conjectured to be $L.end = 71 /Zq. This expression is a function of 
j3 since Zq in general depends on the coupling constant. 

We have run our Monte Carlo simulations at imaginary values of the theta parameter, 6l = 
+i$L (#l £ K.) i n order to avoid the sign problem. Let us briefly introduce the new cluster algorithm 
for imaginary Qi. 

The first part of an updating step with the Wolff algorithm [|l^] for the standard 0(3) sigma 
model without a 6 term consists in choosing a random unit vector r in such a way that every 



dynamical field can be split in a component parallel to r and the rest, <j>(x) = ( 0(x) -rj r + <j>±(x), 

where 0j_(x) denotes the part of <f(x) orthogonal to f. Then the signs of (<j>(x) -rj for all x are 
updated a la Swendsen-Wang as in the Ising model. 



By introducing the above separation for <j>(x) in the expression ( |2.2| ) we can rewrite it as 

Ql(x) = — { (j>(x) -rj (d h2 + d-i - 2 + d 2 -i + d_ 2 ,i ) + 

<j)(x+l)-r) (d 0) -2-do/z) + (<j>{x-l)-r) (d ,2 - <h 2) + 



$(x+2) -r) (do,i -db,-l) + (?(*-2) •?) (^-1-^0,1) } , (2-4) 

where x± 1 means the site at the position one step forward (backward) in the direction "1" starting 
from site x and the notation dtj stands for the 3x3 determinant (the three components for each 
vector are unfold along the rows) 



dij = det I (f(x+i ) ) . (2.5) 
In this fashion the theory at each updating step looks like an Ising model in the presence of an 



external local magnetic field h(x) because the 6 term in Eq.([2~4|) is linear in yp-rj. One can 
readily derive that the magnetic field at site x is 

h(x) = -Jfa${ x ) ' 7 \ {d\,2 + d-i-2 + d2-i+d-2,\+d-i-i-2 + d-i + 2-i + 

d\\+2 +^1-2,1 +^2,2-1 +^2+1,2 + ^-2,-2+1 +^-2-1,-2 ) • (2-6) 



(and analogous terms in (2.6)) are the straightforward generalization of the above defini- 



tion (2.5) when the site is obtained by shifting two steps ( i plus k ) from the original position x. 
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Figure 2: Behaviour of the triplet mass (in units of the lattice spacing) as a function of 2 . Circles (j8 = 1 .5) 
and triangles (j3 = 1.6) are the data from the simulation at imaginary (0 2 < 0). Each continuous line is 
the result of the analytic continuation described in the text and the dashed lines enclose the boundary of its 
error. 



We have converted the original theory in an Ising model in the presence of a local external 
magnetic field (which changes at each updating step and therefore must be recalculated at each 
step). In the literature there are two algorithms expressly introduced to update Ising models in the 
bosom of magnetic fields: the Lauwers-Rittenberg [ |l9| ] and the Wang [ p0| ] methods. We tested 
the performances of both algorithms, compared their decorrelation times and decided for the latter. 
The cluster construction was tackled by the Hoshen-Kopelman algorithm [21]. 



3. Analytic continuation and determination of en d 

Operators which couple the vacuum with the singlet or triplet states can contain an arbitrary 
number of fundamental fields since the model is not parity invariant for 0/0. We only present the 
results for the triplet particle. The following operators were used 

<?i(x)=#(jc), ~$ 2 (x) = -i(f(x)x(f(x+l). (3.1) 

Then we calculated the related wall operators by averaging over the x\ coordinate (as usual L is the 
lattice size), ^ ; (x 2 ) = j£ X] ~&i{x) for i = 1,2. 

To single out the correct parity mixture for the physical particle and to clean the signal from 
possible excited states, we extracted the triplet mass mj by using the variational method of Ref . [ |22| ] 
where the mass is obtained from the largest eigenvalue of the correlation matrix {Wj(x2)Wj(0)) — 
(W i) (W j) . In Fig. ^| the results for the triplet mass are shown for two values of jS. The analytic 
continuations in this figure were done by using a ratio of degree 2 polynomials; analogous results 
are obtained from other functional forms (we checked this statement by using degree 4 and 6 poly- 
nomials). The endpoint (where the mass vanishes) equals (O^end) 2 = (0end/Zg(j6)) 2 where enc j is 
conjectured to be equal to %. On the other hand the renormalization constant Zq can be evaluated 
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Figure 3: Data for (gz.) start at +1 at the 0-th Heat-Bath step and then they go down until reaching a 
plateau. The horizontal line and grey band are the value and error respectively of Zg(/3 = 1.5). 



with the non-perturbative numerical method of Ref. [17]. An example of such an evaluation is 
shown in Fig. ||. Summed up briefly: a classical instanton (with topological charge +1) is put by 
hand on the lattice and then it is heated by applying 100 updating steps (we used Heat-Bath steps 
on the conventional 0(3) non-linear sigma model without a 6 term since Zq cannot depend on 
6). The value of Qi is measured while the continuum topological charge Q is continuously moni- 
tored (by 6 cooling steps after every Heat-Bath updating) to be sure that the background charge is 
not changed. This procedure was repeated 10 4 times for both values of j8. The average of <2l on 
configurations that lie in the topological sector + 1 yields Zq. 

Table 1 



p 


L 


(0L,end) 2 




$end 


1.5 


120 


111(5) 


0.299(18) 


3.15(20) 


1.6 


180 


94(5) 


0.313(12) 


3.03(14) 



Then, once we know Zq and ((^end) 2 the value of d end can be extracted. The results are shown 
in Table 1 . The lattice sizes were chosen large enough to guarantee the absence of relevant finite 
size effects (we imposed L/E, = L • amj > 10). 2 • 10 5 propagators were measured on as many 
independent configurations for all values of j8 and 6l- The results of the last column are in fair 
agreement with the conjecture that predicts enc j = %. By averaging over the results for both /3 
(assuming gaussian errors) we obtain that the endpoint for the triplet mass is en d = 3.07(11). 



4. Conclusions and outlook 



We are studying the spectrum of the 2D 0(3) non-linear sigma model with a 6 term as a 
function of this parameter. The model contains a triplet and a singlet states whose mass depend 
on 6 as shown in the sketch of Fig. [l[ In the present progress report we have given clear evidence 
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that the triplet mass indeed vanishes when 6 becomes %. To do so, we have simulated the model 
at imaginary 6 and then extrapolated the results to real 8. The extrapolation always indicate that 
the mass tends to vanish at an endpoint. Our calculation yields en d = 3.07(11) in good agreement 
with the conjectured prediction en d = 

A new fast cluster algorithm was purposely introduced to simulate the theory at imaginary 6. 

We are planning to improve the statistics by studying the model at other values of j8 and 
other lattice regularizations of the topological charge as well as by increasing the precision in the 



evaluation of Zq which is the largest source of error in the calculation of d en d (see Ref. [14]). We 
will also extend the analysis to the singlet mass and the determination of C . 

We emphasize the good performance of the analytic continuation in our study. It is important 
for that to have got Monte Carlo data within a wide range of values of 6l, (#l = —iOl £ [0, 10]). 
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